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Abstract 

We discuss striped phases as a state of matter intermediate between two extreme 
states: a crystalline state and a segregated state. We argue that this state is very 
sensitive to weak interactions, compared to those stabilizing a crystalline state, and 
to anisotropies. Moreover, under suitable conditions a 2D system in a striped phase 
decouples into (quasi) ID chains. These observations are based on results of our studies 
of an extension of a microscopic quantum model of crystallization, proposed originally 
by Kennedy and Lieb. 

1 Introduction 

Striped phases, specific phases having quasi-one- dimensional structure, are ubiquitous. They 
have been observed in numerous experiments carried out in a variety of systems, among which 
there are physisorbed monolayers on metallic surfaces [I] or ultrathin magnetic garnet films 
[21 [3] . They have also been found in theoretical analyses, carried out by analytic or numerical 
methods, of various quantum and classical model systems. In the class of quantum models, 
one can find results pertaining to the Hubbard model [UEIE], t-J model [HE], XY model [9|, 
spinless fermion model [TUl ITT] , spinless Falicov-Kimball model [T2l [T3] . It has been argued 
in a number of works that some essential physics of systems displaying striped phases can 
be grasped by a kind of coarse-grained descriptions that result in classical lattice-gas models 
(Ising-like models or 0(n) models) with competing two-body interactions. Typically, in 
these models short-range attractive (ferromagnetic) interactions compete with long-range 
repulsive ones [HI [151 EH1 EH EH HH]- Interestingly, despite the difficulty of determining 
ground states of such systems, in recent years a few rigorous results have been published 

[23121]. 

Perhaps the most notable striped phases are those observed in doped layered perovskites, 
which under suitable conditions become high-temperature superconductors [221 123]- The 
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nature of striped phases in these materials, in particular their competition with a supercon- 
ducting state, is still vigorously debated. 

The purpose of this paper is threefold. Firstly, we argue that the problem of the formation 
of striped phases falls naturally into the context of an extended crystallization problem, that 
is not only the existence of crystalline phases is addressed but also the process of their 
deterioration into a segregated state (described below) as a result of weakening of the forces 
stabilizing a crystalline state. A crystalline phase and a segregated one might be thought of 
as two extreme states of matter. Presumably, striped phases appear in this process as the 
last but one stage of a sequence of phase transformations that starts with a crystalline state 
and ends with a segregated one. 

Secondly, we point out that the stability of striped phases might be rather fragile. Very 
weak interactions, in the hierarchy of the interactions present in the system under consid- 
eration, can change the character of the striped pattern or destabilize completely striped 
phases. Moreover, even a weak anisotropy of the system chooses the direction of stripes in 
stable striped phases. 

Thirdly, the appearance of striped phases might signal a dramatic destruction of corre- 
lations in direction of stripes. An anisotropy of the system favors this effect. Specifically, 
we show that when our system is in an axial-stripe phase, zero-temperature electron corre- 
lations (given by ground-state one-body reduced density matrix) along the stripes vanish. 
We make an attempt to estimate the range of temperatures in which this property holds 
approximately. Consequently, before a 2D crystal deteriorates completely (to a state of seg- 
regation), it transforms effectively into a quasi one-dimensional structure, destroying some 
of existing 2D long-range orders. 

Our main conclusions and suggestions are based on a rigorous analysis of an extension of 
a microscopic quantum model of crystallization, which was proposed twenty two years ago 
by Kennedy and Lieb [24|. Some additional observations are based on numerical calcula- 
tions. In brief, the Kennedy and Lieb model of crystallization is a two-component fermionic 
system on a lattice that consists of heavy immobile ions and light hopping electrons, and the 
sole interaction in this system is an on-site electron-ion interaction representing a screened 
Coulomb interaction. 

The paper is organized as follows. In the next section we introduce our model of crys- 
tallization. Then, in Section 3, we present some of the ground-state phase diagrams of our 
system, discuss possible scenarios of deterioration of a checkerboard-like 2D crystal, empha- 
sizing the role of striped phases and the effect of an anisotropy of electron hopping. After 
that, in Section 4, we focus our attention on axial-stripe phases, specifically on electron cor- 
relations in such phases. Finally, in Section 5, we summarize our results. Spectral properties 
of axial-stripe phases are presented in Appendix. 

2 The model of crystallization 

The model of crystallization proposed by Kennedy and Lieb [21] is composed of two fermionic 
subsystems: light hopping electrons and heavy immobile ions. The electrons are represented 
by spinless fermions (spin does not play any role in our considerations), described by creation 
and annihilation operators of an electron at a state localized at site x of the underlying 
lattice, c+, c x , respectively, satisfying the canonical anticommutation relations. The ions are 
described by collections of pseudo-spins {s x } x£A , called the ion configurations; s x = 1 if the 
site x is occupied by an ion and s x = — 1 if it is empty. The pseudo-spins commute with the 
creation and annihilation operators of electrons. 
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There is neither a direct interaction between mobile electrons nor between immobile 
ions. The electrons energy is due to hopping (typically a nearest-neighbor (n.n.) hopping 
is assumed) with t being the n.n. hopping intensity (without any loss of generality we can 
fix the sign of t, t > 0), and due to a screened Coulomb interaction with the immobile 
ions, whose strength is controlled by the coupling constant U > (concerning the sign see a 
comment in the text below). The Hamiltonian of the system reads: 




where (x, y)\ means that the sites x, y constitute a pair of n.n. sites. In the sequel, we shall 
limit our considerations only to the case in which the underlying lattice is a square lattice. 

It is worth to emphasize here that although the ions are not moving due to the dynamics 
of the system, their configurations are not frozen or random. On calculating the canoni- 
cal partition function, one takes a trace over electronic degrees of freedom and sums over 
all configurations of the ions. It is this sampling of possible configurations that produces 
correlations between ions, which in turn may lead to a crystalline arrangement of the ions. 

The Hamiltonian H fk is widely known as the Hamiltonian of the spinless Falicov-Kimball 
model, a simplified version of the Hamiltonian put forward for describing electronic subsys- 
tems of some solids in [25]. To the best of our knowledge, the spinless-fermion Falicov- 
Kimball model is the unique system of interacting fermions, for which the existence of long- 
range orders (periodic phases) and phase separations have been proved. The majority of 
rigorous results refers to the ground states of the model and holds only in the strong-coupling 
regime, i.e. for sufficiently small \t/U\, and for particle densities satisfying specific condi- 
tions (some more details and comments can be found below). Interestingly, the first proofs 
of the existence of a periodic phase, which in 2D is the so called checkerboard phase with 
the electron and ion densities equal 1/2, hold for any U: at zero temperature (26J, [23] and for 
sufficiently low temperatures [24] . It is the checkerboard phase that in our considerations 
is identified with an initial crystal, which is given a possibility to deteriorate. A review of 
results pertaining to ground-state or low-temperature phase diagrams and an extensive list 
of relevant references can be found in [271 [28] ) . 

According to the state of art, a quite general analysis of ground-state phase diagrams, 
that made possible proving the existence of numerous periodic phases and mixtures (phase 
separated states) of such phases, is feasible only in the mentioned above strong-coupling 
regime, and when the densities of the ions and electrons, pi and p e , respectively, satisfy a 
specific condition. Either it is the condition of neutrality, pi = p e , if the electron-ion inter- 
action is attractive, or it is the condition of half-filling, pi + p e — 1, if that interaction is 
repulsive. The two cases are related by a unitary transformation: a hole-particle transfor- 
mation (for definiteness, we assume in our considerations that the electron-ion interaction is 
repulsive and the system is half-filled) . In the specified above regime, it is possible to derive 
explicitly, in the form of a convergent power series with respect to the small parameter \t/U\, 
an effective interaction between the ions (29J, [30J [311 [32] . The components of this expansion 
constitute many-body finite-range lattice-gas interactions, where the number of interacting 
bodies and the range of interaction grow without bound with the order of those components. 
Then, it is possible to construct rigorously the ground-state phase diagram of the effective 
interaction truncated at certain order, and after that derive information on the ground- 
state phase diagram of the complete quantum system under consideration [29j [301 EU 132] . 
The ground-state results obtained in the described way can be extended to sufficiently low 
temperatures [52] . 
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The emerging phase diagram is very reach; besides a few periodic phases J2H ESI EH ES] 
and various phase-separated states [331 El], that can be determined rigorously, it con- 
tains most probably infinitely many periodic phases and mixtures of periodic phases (see 
[36| E3 EH EH] for the results obtained by means of the method of restricted phase dia- 
grams). However at half- filling and large U, it does not contain any segregated phases, that 
is thermodynamic mixtures of the completely filled with ions phase (pi = 1) and the ion void 
phase (pi = 0) with suitable electron densities. This is in agreement with rigorous results 
of [101 HU H2] , which state that a segregated phase is stable for any total density away from 
half-filling if U is large enough. Approximate results of [361 EH EH] suggest that whatever 
U > 0, there is no segregated phases at half-filling. 

We insist, however, that a crystalline state, like the half-filled checkerboard phase, should 
have a possibility to deteriorate into the half-filled segregated phase, considered as a final 
stage of a deterioration of a crystal, for any U. The standard Falicov-Kimball model cannot 
encompass such a "process": there is neither the stable half-filled segregated phase nor a 
parameter to drive such a process. The reason for which half-filled segregated phases are 
missing in the standard spinless Falicov-Kimball model, in the strong-coupling regime, might 
be attributed to the fact that the leading interaction term (whose order is 0(\t/U\ 2 )) of the 
effective interaction expansion is a n.n. repulsive lattice gas (or in terms of pseudo-spins 
- an Ising antiferromagnet), irrespectively of the sign of U. In order to stabilize the half- 
filled segregated phase and to provide the standard model with a suitable control parameter, 
we extend it by adding a small (second order as compared to the electron-ion interaction) 
short-range attractive interaction between the ions. For simplicity, we choose its range to be 
limited to the distance between next-nearest-neighbor (n.n.n.) sites. The physical source of 
this interaction might be identified as van der Waals forces, which despite their weakness are 
known to play an important role in stability of crystals. Therefore, we extend the standard 
spinless Falicov-Kimball model by adding to the Hamiltonian H FK the term H vdW , 



H v dw = — SxS v ~ Yq ^ SxSy, ' ' 

with 

W = -2t 2 + tV e = t 4 e, (3) 

where (x,y)2 stands for a pair of n.n.n. sites. The parameters to and e, which vary the 
strength of n.n. and n.n.n. couplings in H v dw-, respectively, are the control parameters of 
the phase diagrams discussed in the sequel. The specific form of W, given in (j3J), and the 
coefficients that stand by the parameters W and e in (J2]) stem from our analysis of the 
ground-state phase diagram of the 2nd order effective interaction (see [13] for details). In 
those diagrams, u = e = is the point of coexistence of the checkerboard phase and the 
half-filled segregated phase. Despite the weakness of the n.n.n. term in H v aw, it will become 
clear in the discussion of phase diagrams that follows, that extending the range of H V dw 
beyond n.n. enables us to make interesting observations pertaining to a transition from a 
crystalline state to a segregated one. 

To take into account important effects of anisotropy of electron hopping, we differentiate 
between hopping in the vertical and horizontal direction by introducing the corresponding 
hopping intensities t v > and th > 0, and the anisotropy parameter 7 such that t v = y/ythi 
with < 7 < 1 [43J (this somewhat unusual definition of the anisotropy parameter is dictated 
by the simplicity of the expression for the effective interaction). Then, in the case of the 
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hole-particle invariant system, the effective interaction (in the units of U) between the ions, 
up to order 4, reads: 



t 2 3t 4 3 4 W 

7 t 4 + 

1 Id 8 



s x s y + 



' t 2 3 , 2 3t 4 W 
7 4~8 7t ~ 7 16 + T 



(x,y)i,v 



■ 3^_£' 
7 1~6 16 



+4 

t 
16 



(4) 



where t = t^/C/, (x, y) n ,h ((%, y)n,v) stands for pairs of nth-order nearest-neighbor sites (n.n.- 
lst order, etc) in horizontal (vertical) directions, and sp denotes a product of four spins 
whose sites constitute an elementary square, P, on the lattice. 

Since we are working with a truncated effective interaction, we have to assign an order 
to the deviation of the anisotropy parameter 7 from the value 1 (the isotropic case). For 
this purpose we define the anisotropy order, a, and the new anisotropy parameter, f3 a : 7 = 
1— j3 a t a , p a > 0. In order to analyse anisotropy effects with the effective interaction truncated 
at certain order, the anisotropy order has to be suitably adjusted. In our work, we analyze 
the effective interaction truncated at 4th order, hence the weakest admissible deviation from 
the isotropic case corresponds to a = 2. 



3 Ground-state phase diagrams 

The ground-state phase diagrams constructed with effective interaction (T4J) are shown in 
Fig. [TJ The essential point is that they can also be thought of as ground-state phase diagrams 
of the full quantum system, provided that the coexistence lines of the phases are interpreted 
as strips whose width is of the order 0(\t/U\ 2 ). Inside these strips, we cannot claim the 
stability of any phase, but we can exclude the stability of some phases. Studies carried out 
by approximate methods suggest that these strips may accommodate infinitely many phases 
[361 S3 EH] • To unveil these phases rigorously, it is necessary to construct phase diagrams 
according to the effective interactions truncated at higher orders, but this task becomes 
quickly hardly feasible. 

Because of the hole-particle invariance, in all the phases of the phase diagrams shown 
in Fig. [T]the densities of the ions and the electrons are 1/2; off the hole-particle symmetry 
case, many other phases of different than 1/2 densities can be proven to be stable [43J. 
Suppose that initially our system is in a crystalline checkerboard phase, <S C 6, Fig. [TJ Then, 
by decreasing the "van der Waals interaction" coupling u, the system is always driven out 
of the crystalline state. However, the terminal state and the way it is attained depend 
strongly on variation of the tiny coupling controlled by e. It is worth to note here that the 
complete n.n.n., 4th order, effective interaction consists of a contribution from Hpx and from 
H V dw] it is attractive if e > 3, irrespectively of the value of the electron hopping anisotropy 
controlled by 02, and then it stabilizes the checkerboard phase. Concerning a terminal state, 
one can observe that if e does not decrease (while u decreases), then the terminal state 
is always the segregated phase, which is, in this case, a fifty-fifty mixture of the uniform 
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Figure 1: Ground-state phase diagrams in the case of the hole-particle invariant system, 
according to the effective interaction The representative ion configurations of the dis- 
played phases are shown in Fig. [2J The superscript v by the symbols of phases means that 
only phases with vertically oriented configurations of ions are stable. S seg stands for the 
segregated phase, which is a mixture of two uniform phases: S + (where s x — 1 at every site 
x) and S- (where s x = —1 at every site x). 
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Figure 2: Representative ion configurations of the phases considered in the paper. S c b and 
Sdi are examples of checkerboard-like crystals, Sdd - of diagonal-stripe phases, <Si, «S 2 , S3 - 
of axial-stripe phases. 1S3 is absent in the phase diagrams of Fig. [TJ 

phases S + and *S_. If, on the other hand, e does not increase, then whether the segregated 
phase is attained or not depends on the rate of its variation as compared with that of u. 
In particular, if e decreases sufficiently fast, then the process of deterioration of the crystal 
terminates in one of the axial-stripe phases <Si, <Sf, S2 , S%. The state of segregation is not 
attained. This is reminiscent of the results of qualitative analyses of classical lattice gases 
with competing interactions, repulsive long-range dipole-dipole or Coulomb interactions and 
attractive short-range interactions [HI \%5[ EESj. 

Now, concerning the sequence of intermediate phases visited by the system before the 
terminal state is attained, we can see in Fig. [1] only a few elements of this sequence (see 
the remark above concerning the coexistence lines). Apparently, two cases can be realized. 
In the first one, the only intermediate phase is a mixture of S c b, S+, and <S_. This occurs, 
for instance, if the value of e in the initial state is large enough and is kept constant. In 
this case, the system being initially in a crystalline state attains the terminal segregated 
phase directly via a 1st order phase transition [171 IS]. in the second case, before a terminal 
state (the segregated phase or an axial-stripe phase) is attained, the system visits numerous 
intermediate phases; at the end of this sequence one always finds axial-stripe phases. This 
case is realized, for instance, if e is sufficiently small and kept constant. Then, the n.n.n., 
4th order, effective interaction becomes repulsive and frustrates the interactions stabiliz- 
ing the checkerboard phase. As a result, as to decreases, first the ions form dimers which 
are distributed in a checkerboard-like manner in the phase Sdi, then the the particles form 
extended strings, i.e. completely filled lattice lines, like in the axial-stripe phases, Si and 
S3. The latter succession of phases is reminiscent of those observed in many other systems 
with competing interactions (see for instance [HI US] for results of computer simulations 
and [21J for rigorous results) and conforms to the following observation: purely repulsive 
interactions tend to proliferate as small as possible aggregates of particles and arrange them 
in checkerboard-like crystalline patterns, while attractive interactions tend to make the ag- 
gregates of particles as large as possible, which results eventually in segregated-like phases 
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with mesoscopic regions completely filled with particles or completely empty. Let us em- 
phasize again that, when a transition from the checkerboard crystal to a segregated phase 
or a striped phase is accomplished, the sequence of visited phases is, most probably, much 
longer. The phases visible in the diagrams of Fig. [TJ being a part of this sequence, point out 
only to some tendencies. Apparently, in a process of deterioration of a crystal, striped phase 
constitute a stage just preceding a transition to a segregated phase. 

Finally, let us consider an impact of electron hopping anisotropy on the phase diagrams. 
For defmiteness, we set the vertical hopping weaker than the horizontal one (i.e. 02 > 0). 
Any hopping anisotropy breaks the rotational symmetry of the system and out of the phases 
that are not invariant with respect to all rotations, select some with a specified orientation. 
For instance, among the axial-stripe phases it stabilizes those stripes that are oriented in the 
direction of a weaker hopping (vertical in our case). This effect is clearly visible in the phase 
diagrams of Fig. [TJ where for a nonzero anisotropy of n.n. hopping only vertically oriented 
striped phases remain stable. Let us emphasize that this property holds not only for a 
truncated effective interaction for which the phase diagrams are actually constructed but for 
the full quantum system as well. In [12], it was demonstrated that an analogous effect occurs 
also for diagonal-stripe phases in the presence of an anisotropy of n.n.n. electron hoppings. 

The effect of orienting axial-stripe phases in the direction of a weaker hopping has been 
found also in studies of a Hubbard model in the framework of real-space Hartree-Fock ap- 
proximation [6]. 

4 Specific properties of axial-stripe phases 

All the conclusions of the previous section are based on a rigorous analysis of ground-state 
phase diagrams of our system, for sufficiently large electron-ion coupling and for a half-filled 
system, with the main result referring to the stability of striped phases. In this section, we 
look into properties of the axial-stripe phases to see how important are the above conditions 
for the drawn conclusions, whether the effect of orienting axial stripes in the direction of 
weaker hopping, seen in the diagrams of previous section, might persist beyond the regimes 
of strong coupling and half-filling. Moreover, it should be interesting to unveil the electronic 
properties of axial-stripe phases, which have been left untouched by the analysis of previous 
section. 

Numerical calculations of the thermodynamic functions considered below are based on 
exact diagonalization of the electron subsystem in axial striped phases Si, S2, and £3 (the 
corresponding spectra and some of their properties are reproduced in Appendix). In this 
section we present only some properties, verified for the phases Si, S 2 , and S 3 (see Fig.[2J the 
lack of subindex "v" indicates that we consider both the horizontally and vertically oriented 
stripes). We think that the presented properties are characteristic for the whole family 
of the axial-stripe phases. To support our conclusions, we display only plots of discussed 
quantities for the phase £3 (see Fig. [2]), which is not present in the displayed here hole- 
particle symmetric phase diagrams (it is visible in diagrams off the hole-particle symmetry 
[T31 H5]). but appears to be a good representative of the family of axial-stripe phases. 

We start with comparing the ground-state internal-energy densities of vertically and hor- 
izontally oriented axial-stripe phases of a half-filled system, for a small and a large electron- 
hopping anisotropy, see Fig. [31 Let us recall that according to the phase diagrams displayed 
in previous section, it is the phase whose stripes are oriented in the direction of weaker 
hopping (vertical in our case) that has a lower energy. This effect occurs in a situation, 
where all the phases in the phase diagrams constructed with the method used in our work 
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(strong-coupling expansion of half-filled system) are insulating. The reason is that for U so 
large that the expansion is convergent, a gap opens at the Fermi level of any stable periodic 
phase [27]. 




Figure 3: Energy density (per particle), e, and the width, 5, of the gap at the Fermi level in 
the half-filled phases and S%, as functions of U, for two values of the anisotropy parameter 

7- 

Now, according to our calculations carried out for axial stripes, the vertically oriented 
stripes seem to have lower energy than those perpendicular to this direction, for all values of 
U > 0. Of course, we cannot exclude numerically that this holds only above some positive 
value of U, but in any case this value is much smaller than the values of U admissible in 
the expansions referred to in Section 2. It follows from dispersion relations (reproduced in 
Appendix) that for sufficiently large U, in both systems, of vertical and horizontal stripes, 
there is a gap at the Fermi level. This gap opens for smaller values of U and is larger for 
vertical stripes (those along the direction of weaker hopping), than the horizontal ones. But 
the relative stability of vertical stripes persists even for the values of U, for which there is no 
gap at the Fermi level. The above observations shine more light on the stability of stripes 
oriented in the direction of weaker hopping, relative to perpendicular stripes: the effect 
found in the strong-coupling phase diagrams, when the discussed phases are insulating, may 
persist also when they are metallic. 

However, if we abandon the condition of half- filling the situation is different. In Fig. 0] 
we display the plots of the ground-state internal-energy densities of vertical and horizontal 
stripes, for a small and a large electron hopping anisotropy, and for electron densities p e = 1/2 
and p e = 1/3. In the cases displayed in Fig. H]a,b,c there is no gap at the Fermi level, whether 
the stripes are vertical or horizontal, while in the case of Fig. HU the gap at the Fermi level 
opens for vertical stripes (because of sufficiently strong hopping anisotropy). Nevertheless, 
in all the cases there is a critical value of U above which it is the horizontally oriented 
phase that has lower energy. That is, for large U but away of half-filling, the vertical stripes 
exchange their stability with horizontal stripes. 

Some of the spatial correlations properties of the electron subsystem are determined 
by off-diagonal matrix elements (in position basis {|r)}) of the one-body reduced density 
operator p^ (hereafter called the correlation functions), which can be expressed by the 
eigenf unctions in position representation. If the Fermi level is located at an upper edge of 
a band, which does not overlap with the bands of higher energy, then in the limit of an 
infinite system, the correlation function depending on two lattice positions r and r ; assumes 
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a) p e =l/2, y=0A9 b) p e =l/2, ^0.01 




Figure 4: Energy density (per particle), e, in the phases S3 and 1S3 for two values of electron 
density off the half-filling, and for two values of the anisotropy parameter 7. Only in case 
d) a gap of width S opens at the Fermi level of <Sf. 



the form 



7T 7T 

(r^V) = J dh J ^X(r,k)X*(r',k)exp(ik(r-r / )), 



-7T — 7T 



where the function K(r, k) is given in terms of the components (in the plane-wave basis) of 
eigenvectors corresponding to completely filled bands. As a matter of fact, in the case of 
axial-stripe phases these components, hence the function K(r, k), do not depend on the wave 
vector component in the direction of stripes (see Appendix). Let us consider for definiteness 
the case of vertically oriented stripes. Then, in the vertical direction, that is for r' = r + mj, 
with m being an integer and j standing for a lattice translation vector in the vertical direction, 
the correlation function can be written as 

7T 

dk h K(r,k h )K*(r',k h ) / dk v exp (imk v ), 



where the last integral vanishes identically. 

Let us consider now half-filled axial-stripe phases for fixed U. A sufficiently strong 
anisotropy (a sufficiently weak hopping along the stripes) opens a gap at the Fermi level, 
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and then the chains perpendicular to the stripes become independent. Therefore, for suffi- 
ciently large anisotropy the ground-state extensive thermodynamic quantities of the whole 
2D system become equal to those of ID chains perpendicular to the stripes (i.e. with zero 
hopping in the direction of stripes). This effect is clearly visible in Fig. [5^i, where for the 
half-filled, vertically oriented, axial-stripe phase S% we show the plots of ground-state inter- 
nal energy densities versus U, for different anisotropies. Starting from sufficiently large U, 




Figure 5: Energy density (per particle), e, and relative energy density, Ae = |(e(7) — 
e(0))/e(0)|, for three nonzero values of the anisotropy parameter 7, in the phase S%. The 
vertical lines are located at the values of U, where the gap opens. The highest (continuous) 
curve in part a) corresponds to the ID case. 

each plot of a 2D system coincides with the plot for the corresponding ID system. To see 
more clearly that the plots merge exactly at the value of U for which the gap at the Fermi 
level opens, we have plotted also the relative differences of internal-energy densities Ae (see 
Fig. [5b and its caption). 

It is quite clear that thermal fluctuations will destroy the ideal independence of chains at 
zero temperature. However, it is interesting to learn how effective they are, by calculating, 
for instance, the Helmholtz free-energy- density relative differences, A/, versus temperature 
T (see Fig. Ek,b and its caption). We display the corresponding plots for three different 
anisotropies and two values of U. Let us note that the larger is U and the stronger is the 
anisotropy the larger is the gap at the Fermi level in stripes oriented in the direction of the 
weaker hopping. For the smaller of the U values (Fig. [6^), the gap opens only for sufficiently 
strong anisotropy. In the isotropic case (the highest curve) there is no gap, the chains 
perpendicular to the stripes are dependent even at T = 0, so the plot of the corresponding 
Af starts above zero. For the larger of the U values (Fig. [6b) there is a gap for any 
anisotropy and all the curves start at zero. It is apparent that for a specified tolerance for 
the deviation Af (of a few per cent of the maximal value), the larger is the gap the higher is 
the temperature above which this tolerance is exceeded. Below this temperature the chains 
can be considered as approximately independent. 



5 Summary 

We have studied an extension of a microscopic quantum model of crystallization, proposed by 
Kennedy and Lieb, in which light hopping electrons interact on-site only with heavy immobile 
ions. The original model has been extended by a weak short-range attractive interaction 
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a) p e =2/3, U=4.0 
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b) p e =2/3, U=5.0 
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Figure 6: Relative difference of free energy densities, Af = |(/(7) — /(0))//(0)|, for three 
values of the anisotropy parameter 7, in the phase S% with two values of U. 

between ions, whose purpose is to mimic the effect of van der Waals forces. In the framework 
of the extended model, the stability of striped phases has been analyzed on a rigorous basis. 
Moreover, we have looked closely into properties of a special class of striped phases, the 
axial-stripe phases, particularly into properties of their electron subsystem under presence 
of a hopping anisotropy. It is known already that an anisotropy of the electron hopping 
has a strong impact on the properties of electron systems |50J ; even a half- filled system of 
free electrons in a constant potential (no gap at the Fermi level) turns into insulator, in 
the direction of a weaker hopping, as soon as a hopping anisotropy is "switched on". This 
is a consequence of an exponential decay of correlations in this direction, for any nonzero 
anisotropy. In the model considered here, the effects of hopping anisotropy are even more 
striking. We have demonstrated that for sufficiently large anisotropy (depending on the 
value of U), the stable axial stripes (i.e. oriented in the direction of a weaker hopping) 
decouple into ID chains, which are perpendicular to the stripes. The reason is that (for a 
given U) a sufficiently strong anisotropy opens a gap at the Fermi level, and this in turn 
leads to vanishing of electron correlations along the stripes. Consequently, while the ion 
subsystem develops a 2D long-range order, the electron subsystem may develop only a ID 
long-range order, and the compound 2D system behaves like a collection of chains. It is 
tempting to suggest that an analogous effect may occur in many other systems developing 
stripes, observed experimentally or studied theoretically, some of which were mentioned 
in the Introduction. The appearance of stripes in some degrees of freedom may signal a 
significant reduction of correlations along the stripes among other degrees of freedom of a 
compound system. 



6 Appendix 

Here we present some spectral quantities for electrons in periodic potentials given by vertical 
ion configurations S%, <S| and Sf, under periodic boundary conditions. We denote the 
dispersion relations as Am (k) and the gaps widths as 8 m , where m labels configurations 
(m = 1,2,3) and I counts the bands and gaps from bottom to top. The corresponding 
quantities for horizontal configurations can be obtained by exchanging t v and th in the 
formulae below. 
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First, for all the considered cases we define e„ as 



Sy — COS (kyj • 



Then, for S% 



A?(k) = e v + u±A f 
5 2 = \U\-4\t v \, 

where Eh = 2t h cos(k h /2), u—U/2 and = \/e\ + u 2 ; 



for <Sf 



A?(k) = e v + u + 2^^cos (^ + 2 fj. 



1 = 1,2,3; 



r(0 



-4^1 + 2^/3(^ + ^)8^^+ ;:; 
/ = 1,2; 



£ , 2vr(/ - 1) 



. u 3 + t 3 h cos(k h ) \ ( u 3 + (-ly-^tH 

where u = U /3, (ph = arccos | — — z^75 — I anc * Vh = arccos 



for &[ 



Ai°(k) = e v + u± J2t 2 + u 2 ±2\t h \A 



(3) 



(3 + ^) 3/2 



[U 2 + 



^ = -4N + + (M + M) 2 - \[t\ + (1**1 - Mf 

5{ 2) = -A\t v \+2^t 2 +u 2 -\t h \ ) j, 



Si 



(i). 



where u — U/2,A h — ^t 2 cos 2 {k h /2)+u 2 . 

It might be helpful to represent the above formulae graphically. We have chosen to plot 
the band edges separated by a gap and the Fermi levels as functions of U. In Figs. [7] and 
[HI these plots are displayed for ion configurations «S| and S% , in all the cases considered in 
previous section. Let {|k)j} be a plane-wave basis of the space of single electron states, la- 
belled by wave- vectors k and index i — 1, . . . , I enumerating the sublattices of the underlying 
periodic ion configuration. Let be the matrix, in the plane-wave basis, of the Hamiltonian 
of a single electron in an external field specified by an axial-stripe configuration of horizontal 
period I (which amounts to the number of sublattices). The components of this matrix are 
labelled by the sublattices of the underlying periodic ion configuration. The matrix has the 
following structure: 

/e« + C! f h2 {k h ) ... f h i{k h )\ 
h,i{k h ) e v + C 2 ... f 2 ,i(k h ) 

\fl,l{ k h) flA k h) ■■■ £v + CiJ 

where the matrix elements fij(kh) are some functions of wave vector component kh, and Cj 
are independent of the wave vector. In this matrix, the only elements which depend on wave 
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a) S3, Y=0.49 
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Figure 7: Band edges and Fermi levels for 
configurations 5f and Solid lines repre- 
sent the band edges that are separated by a 
gap (note that in cases a), b) and d) there is 
no gap between the first and the second band) 
while dashed, dash-dotted and dotted lines - 
the Fermi levels. When a line of a Fermi level 
joins an intersection of two solid lines, then 
it should be continued along the lower of the 
solid lines. These plots correspond to the cases 
shown in Fig. [3] and Fig. 01 





Figure 8: Band edges and Fermi levels for con- 
figuration c>3 with electron density p e = 2/3,. 
Solid lines represent the band edges (note that 
in cases c) and d) there is no gap between the 
first and the second band) and the dashed line 
- the Fermi level. When a line of a Fermi level 
joins an intersection of two solid lines, then it 
should be continued along the lower of the solid 
lines. In the case a) the Fermi-level line coin- 
cides with the upper edge of the second band. 
These plots correspond to the cases shown in 
Fig.EJ 
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vector component k v are the diagonal elements, (h^u = e v + Cj. Then, the solutions of the 
characteristic equation assume the form: A(k) = e v +f(kh), with / being some function of the 
horizontal component of k. Consequently, the equations for eigenvectors a with components 
a ii Si(^k — A(k))ijQ!j = 0, and so their solutions, do not depend on the vertical component 
ofk. 
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